EXTENDED PDF FORMAT 
SPONSORED BY 



www. rndsystcms.com 



A Reconciled Estimate of Ice-Sheet Mass Balance 

Andrew Shepherd et al. 

Science 338 , 1183 (2012); 

DOI: 10.1 126/science. 1228 102 


HIaaas 


This copy is for your personal, non-commercial use only. 


If you wish to distribute this article to others, you can order high-quality copies for your 
colleagues, clients, or customers by clicking here. 

Permission to republish or repurpose articles or portions of articles can be obtained by 
following the guidelines here. 

The following resources related to this article are available online at 
www.sciencemag.org (this information is current as of May 3, 2013 ): 

A correction has been published for this article at: 
http://www.sciencemag.Org/content/338/6114/1539.1.full.html 

Updated information and services, including high-resolution figures, can be found in the online 
version of this article at: 

http://www.sciencemag.org/content/338/61 11/1 183.full.html 

Supporting Online Material can be found at: 

http://www.sciencemag.Org/content/suppl/2012/1 1/28/338.61 11.11 83.DC1 .html 

A list of selected additional articles on the Science Web sites related to this article can be 
found at: 

http://www.sciencemag.org/content/338/61 1 1/1 183.full.html#related 

This article cites 156 articles, 32 of which can be accessed free: 
http://www.sciencemag.org/content/338/61 1 1/1 183.full.html#ref-list-1 

This article has been cited by 3 articles hosted by HighWire Press; see: 
http://www.sciencemag.org/content/338/61 11/11 83.full.html#related-urls 

This article appears in the following subject collections: 

Atmospheric Science 

http://www.sciencemag.org/cgi/collection/atmos 


Science (print ISSN 0036-8075; online ISSN 1095-9203) is published weekly, except the last week in December, by the 
American Association for the Advancement of Science, 1200 New York Avenue NW, Washington, DC 20005. Copyright 
2012 by the American Association for the Advancement of Science; all rights reserved. The title Science is a 
registered trademark of AAAS. 


Downloaded from www.sciencemag.org on May 3, 2013 




RESEARCH ARTICLES 


the assembly. However, in practice, low yields 
were already observed for larger designs (up to 
24,576 nucleotides attempted thus far). Solving 
this challenge may require improvements in struc- 
ture and sequence design, enzymatic synthesis for 
higher-quality strands, optimized thermal or iso- 
thennal (44) annealing conditions, and a detailed 
understanding and perhaps explicit engineering 
of the kinetic assembly pathways ( 8 , 14, 44) of 
DNA brick structures. 

The DNA brick structure, with its modular ar- 
chitecture, sophisticated geometry control, and 
synthetic nature, will further expand the range of 
applications and challenges that nucleic acid nano- 
technology has already started to address — for 
example, to arrange technologically relevant guest 
molecules into functional devices (6, 25, 32-34), 
to serve as programmable molecular probes and 
instruments for biological studies (33, 34, 36), 
to render spatial control for biosynthesis of use- 
ful products (25), to function as smart chug deliv- 
ery particles (37), and to enable high-throughput 
nanofabrication of complex inorganic materials 
for electronics or photonics applications (6, 32). 
The modularity of the brick structure may facil- 
itate rapid prototyping of diverse functional nano- 
devices. Its sophisticated and refined geometrical 
control may enable applications that require high- 
precision arrangements of guest molecules. Be- 
cause the brick structure is composed entirely of 
short synthetic strands (no biologically derived 
scaffold), it is conceivable to make bricks by using 
synthetic infonnational polymers other than the 
natural form of DNA. Such polymers may in- 
clude L-DNA (26), DNA with chemically modi- 
fied backbones or artificial bases, or chemically 
synthesized or in vitro (or even in vivo) transcribed 
RNA. This material diversity may potentially 
produce nanostructures with not only prescribed 
shapes but also designer chemical (or bio- 
chemical) properties (such as nuclease resistance 
or reduced immunogenicity) that would be useful 
for diverse applications requiring the structure to 
function robustly in complex environments, 
such as in living cells or organisms. 
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We combined an ensemble of satellite altimetry, interferometry, and gravimetry data sets using 
common geographical regions, time intervals, and models of surface mass balance and 
glacial isostatic adjustment to estimate the mass balance of Earth's polar ice sheets. We find that 
there is good agreement between different satellite methods — especially in Greenland and 
West Antarctica — and that combining satellite data sets leads to greater certainty. Between 1992 
and 2011, the ice sheets of Greenland, East Antarctica, West Antarctica, and the Antarctic 
Peninsula changed in mass by -142 ± 49, +14 ± 43, -65 + 26, and -20 ± 14 gigatonnes year -1 , 
respectively. Since 1992, the polar ice sheets have contributed, on average, 0.59 ± 0.20 millimeter 
year 1 to the rate of global sea-level rise. 


F luctuations in the mass of the polar ice 
sheets are of considerable societal impor- 
tance, because they affect global sea lev- 
els (1, 2) and oceanic conditions. They occur as 


a consequence of their internal dynamics and 
changes in atmospheric and oceanic conditions 
(3-5). Analysis of the geological record sug- 
gests that past climatic changes have precipitated 
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sustained ice-sheet contributions, in excess of 
10 mm year -1 over millennial time periods (6), 
and the prospect of such changes in the future are 
of greatest concern. Even the modest rises in 
ocean temperature that are predicted over the 
coming century (7) could trigger substantial ice- 
sheet mass loss through enhanced melting of ice 
shelves (8-10) and outlet glaciers (11, 12). How- 
ever, these processes were not incorporated into 
the ice-sheet models that informed the current glob- 
al climate projections (13). Until this is achieved, 
observations of ice-sheet mass imbalance remain 
essential in determining their contribution to sea 
level. 

Satellite geodesy has revolutionized the manner 
in which ice-sheet mass balance is estimated 
(14, 15). Since 1998, there have been at least 29 
ice-sheet mass balance estimates, based variously 
on the satellite techniques of altimetry, interfer- 
ometry, and gravimetry (16). These estimates, and 
their respective uncertainties, allow for a com- 
bined Greenland and Antarctic ice-sheet mass im- 
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balance of between -676 and +69 gigatonnes 
(Gt) year -1 , equivalent to a mean global sea-level 
contribution in the range of +1 .9 to -0.2 mm year -1 . 
However, much of this spread, which is large in 
comparison to other ice-sheet imbalance assess- 
ments (1, 2) and to the estimated rate of global 
sea-level rise (1 7), is due to the brevity of many 
satellite surveys (4.5 years, on average) relative 
to the rate at which ice-sheet mass fluctuates 
(5, 18-20). Because the various satellite methods 
differ in their strengths and weaknesses (14, 15, 21), 
careful consideration ought to make them com- 
plementary. Here, we compare and combine es- 
timates of ice-sheet mass balance derived from 
all three satellite geodetic teclmiques, using com- 
mon spatial and temporal domains, to investigate 
the extent to which the approaches concur and to 
produce a reconciled estimate of ice-sheet mass 
balance. 

Data and Methods 

In this assessment, we use 19 years of satellite 
radar altimeter (RA) data, 5 years of satellite laser 
altimeter (LA) data, 19 years of satellite radar 
interferometer data, 8 years of satellite gravime- 
try data, 32 years of surface mass balance (SMB) 
model simulations, and estimates from several 
glacial isostatic adjustment models, to produce a 
reconciled estimate of ice-sheet mass balance. The 
satellite data sets were developed by using in- 
dependent methods and, in the case of the LA, 
gravimeter, and SMB data sets, through contribu- 
tions from numerous research groups. To enable a 
direct comparison, we reprocessed the geodetic 
data sets with use of common time intervals and 
common definitions of the East Antarctic, West 
Antarctic, Antarctic Peninsula, and Greenland ice- 
sheet (EAIS, WAIS, APIS, and GrIS, respectively) 
boundaries (16). The maximum temporal extent 
of the satellite data sets spans the period 1992 to 
2011, and results from all geodetic techniques are 
available between January 2003 and December 
2008. Unless stated otherwise, all results are 
presented with 1 -sigma uncertainty estimates. 

Ice-sheet surface mass balance. SMB in- 
cludes solid and liquid precipitation, surface sub- 
limation, drifting snow transport, erosion and 
sublimation, and meltwater formation, refreez- 
ing, retention, and runoff. Our estimates of the 
Antarctic Ice Sheet (AIS) and the GrIS SMB are 
derived from reconstructions of the RACM02 
regional atmospheric climate model (22) over 
the period 1979 to 2010, with horizontal reso- 
lutions of 27 (MS) and 1 1 (GrIS) km. RACM02 
has a multilayer snowpack with drifting snow and 
snow albedo schemes (23) and has been evalu- 
ated against in situ temperature, wind, and surface- 
energy balance observations from weather stations 
(24-26) as well as satellite-derived estimates of melt 
extent, mass changes, and drifting snow (5, 26, 27). 
The spatial uncertainty of the RACM02 mean 
SMB has been assessed through comparison with 
310 (GrIS) and 1850 (AIS) in situ observations 
(28). However, temporal fluctuations in snow 
accumulation are poorly resolved in observation- 


al data sets, so we assess the temporal uncertainty 
through comparison with global atmospheric 
reanalyses (16). We also use RACM02 to drive 
a model of MS fim densification (29) for the 
purpose of converting satellite LA observations 
into changes in ice-sheet mass. 

Glacial isostatic adjustment (G1A). GLA of the 
solid Earth is an important contributor to the 
signals observed by satellite gravimetry and, to a 
lesser extent, satellite altimetry (30). The GIA 
must therefore be considered when estimating ice- 
sheet mass balance with either technique. In 
Antarctica, the use of GIA models has in prac- 
tice introduced considerable uncertainty (up to 
1 30 Gt year -1 ) into ice-sheet mass balance estimates 
derived from satellite gravimetry (31-33). There 
are a number of contributory factors to this un- 
certainty, including the scarcity of constraints on 
the evolution of the ice sheet since the Last 
Glacial Maximum (LGM), limited knowledge of 
Earth mechanical properties, and the scarcity of 
near-field relative sea-level and vertical crustal 
motion data with which to evaluate model per- 
formance (34, 35). 

Here, we consider variants of six GIA mod- 
els, and we assess their impact on geodetic ice- 
sheet mass balance estimates. For Greenland, 
where the signal of GIA is relatively small and 
well constrained, we use the Simpson (36), ICE- 
5G (37), and ANU (35) models. For Antarctica, 
we compare the ICE-5G model (39) with two 
recent Antarctic GIA models: the W12a model 
(35, 40) and a version (IJ05 R2) of the H05 model 
(41) updated for this study (16). Both regional 
GIA models incorporate recently improved con- 
straints on the ice-loading history (42—45) that 
suggest that the MS was thinner at the LGM than 
previously thought, leading to a lowering of es- 
timated ice-sheet mass losses since that time 
(40, 45). Although a consequence of this re- 
vision is a potential discrepancy between far-field 
sea-level records and commonly accepted North- 
ern Hemisphere deglaciation models, both of the 
new regional GIA models perform well when 
compared with Antarctic Global Positioning Sys- 
tem (GPS) observations (34), and we conclude 
that these latest solutions are best suited for esti- 
mating MS mass balance. 

Radar and laser altimetry. RA and LA pro- 
vide ice-sheet mass balance estimates through 
measurements of ice-sheet volume change. The 
technique has been applied to both Greenland 
(46-48) and Antarctica (4, 47, 49) and is unique 
in spatially resolving the detailed pattern of mass 
imbalance, with monthly temporal sampling. RA 
provides the longest continuous record of all geo- 
detic techniques (50). Mtimeter measurements of 
elevation change are precise, because they require 
only modest adjustments to account for sensor 
drift, changes in the satellite attitude, atmospheric 
attenuation, and movements of Earth’s surface. 
By far the greatest uncertainty lies in the con- 
version from volume to mass change. In the case 
of LA, this conversion has been performed by 
using an external model of fluctuations in the fim- 
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layer thickness (29, 48, 51). In the case of RA, the 
conversion to mass has been performed by using 
a prescribed density model and by allowing for 
temporal fluctuations in snowfall in the uncer- 
tainty (52). 

We used European Remote-Sensing (ERS-1 
and ERS-2) satellite and Envisat 35-day repeat 
satellite RA observations to determine changes in 
the mass of the EAIS and WAIS between May 
1992 and September 20 10 (16). Time series of sur- 
face elevation change were developed at 39,375 
crossing points of the satellite orbit ground tracks 
by using dual-cycle crossovers (49, 53). In total, 
46.5 million measurements were included in this 
analysis, encompassing 74 and 76% of the EAIS 
and WAIS, respectively. The satellites were cross- 
calibrated by considering differences between ele- 
vation changes occurring during periods of mission 
overlap. Elevation data were corrected for the lag 
of the leading-edge tracker and for variations 
in dry atmospheric mass, water vapor, the iono- 
sphere, solid Earth tides, and surface scattering 
(50). The IJ05 R2 model was used to correct 
for elevation changes associated with GIA. Mass 
changes were calculated by using a surface-density 
model with a nominal density for fim (400 kg nT 3 ) 
applied to all regions other than those in which 
changes are assumed to occur at the density of ice 
(900 kg in 3 ) (52). Rates of mass change were 
computed in regions of interest by interpolating 
measurements derived at satellite-orbit crossing 
points and by extrapolating these results to un- 
observed area. To estimate the uncertainty of mass 
trends, we treated the estimated variability of snow- 
fall (49) and the elevation trend variability as 
equivalent sources of uncertainty. This approach 
is used because it has not yet proved possible to 
separate, in the observed elevation change, an- 
nual cycles due to density fluctuations from res- 
idual variations due to signal penetration into 
the fim. 

We used ICESat (Ice, Cloud, and Land Ele- 
vation Satellite) LA observations acquired be- 
tween September 2003 and November 2008 (the 
period of optimal instrument calibration) to es- 
timate changes in the mass of the AIS and GrIS 
(16). AIS and GrIS elevation rates were com- 
puted by four and two different groups, respec- 
tively, using methods that compare surface heights 
measured along repeated ground tracks. This ap- 
proach provides fine along-track resolution with 
high precision (54). However, the ground tracks 
are widely separated at lower latitudes, and the 
elevation data are sparsely sampled in time be- 
cause of the episodic nature of mission cam- 
paigns and the presence of clouds. The elevation 
data were corrected for the effects of GIA by 
using the W12a model in Antarctica and a com- 
bination of models in Greenland, and three groups 
corrected AIS measurements for estimates of the 
systematic bias between mission campaigns (32). 
A variety of approaches were used to isolate ob- 
servations affected by clouds and to interpolate 
elevation rates between ground tracks. Elevation 
rates were adjusted for the effects of short-term 


fluctuations in fim thickness by using models 
driven by either regional climate model predic- 
tions (29) or by remotely sensed estimates of 
temperature (48, 51). The error budget was cal- 
culated from uncertainties in the corrected height 
measurements, in the correction for change in 
fim thickness, and in the estimated SMB. The 
mass-change estimates reported here are the arith- 
metic average of those obtained by the different 
groups. 

Input-output method (IOM). The 10M quan- 
tifies the difference between glacier mass gained 
through snowfall and lost by sublimation and 
meltwater runoff and the discharge of ice into the 
ocean. The approach has the advantage of al- 
lowing changes in SMB and ice dynamics to be 
examined separately at the scale of individual 
glacier drainage basins (5) and has been used in 
numerous assessments of AIS and GrIS mass ba- 
lance (18, 55-57). Although earlier IOM studies 
used representations of SMB developed from 
guided interpolation of sparse ground observations 
(58-60), regional atmospheric climate models 
(5, 61) are now used because they provide sub- 
daily predictions at high spatial resolution that 
are independent of the in situ observations. When 
evaluated against such data, SMB model errors 
are found to range between 5 and 20%, depend- 
ing on the basin size and location, with propor- 
tionately the largest uncertainties occurring in 
regions of extreme (low or high) precipitation, 
strong melting, or where the model resolution is 
too coarse. Quantifying ice-stream discharge re- 
quires measurements of ice velocity and thickness 
at tire grounding line. Ice-sheet velocity snapshots 
have been widely measured by using interfero- 
metric synthetic aperture radar (InSAR) with high 
(<3%) accuracy (62, 63) and relatively low (an- 
nual or longer) frequency. The thickness of many 
ice streams has been directly measured by using 
airborne radar with high (~10 m) accuracy (64). 
Nonetheless, there are many ice-sheet outlet gla- 
ciers for which such data do not exist; in these 
regions, less accurate methods are used to cal- 
culate thickness with uncertainties in the range of 
80 to 120 m (18, 55). Lastly, where thinning rates 
are large, the temporal evolution of ice thick- 
ness should be accounted for (65). 

We used the IOM to determine mass changes 
of AIS and GrIS drainage basins between 
January 1992 and June 2010 (16). These results 
were derived according to the method of (5 7) and 
updated to include more recent data sets. SMB 
estimates for Greenland as reported in (66) were 
extended to the end of 2010. For Antarctica, the 
SMB estimates of (61) were used, with an up- 
dated uncertainty estimate (28). Ice discharge 
rates were updated by using new ice-thickness 
measurements in the Bellingshausen Sea sec- 
tor, Wilkes Land, and the Amundsen Sea sec- 
tor; at the grounding line of Filchner Ice Shelf; 
and at Byrd and Lambert glaciers. Direct mea- 
surements of ice thickness are now available for 
nearly all WAIS ice streams. The IOM inven- 
tory, including measured and derived thick- 


nesses, now encompasses 64, 79, 96, and 93% of 
the APIS, EAIS, WAIS, and GrIS, respectively; 
the remainder is assumed to have no loss due to 
ice dynamics. 

Gravimetry. The Gravity Recovery and Climate 
Experiment (GRACE) satellite mission has 
allowed fluctuations in ice-sheet mass to be esti- 
mated through measurement of their changing 
gravitational attraction (32, 67-69). Advantages 
of the GRACE method are that it provides 
regional averages without the need for interpola- 
tion, measures the effect of mass fluctuations 
directly, and pennits monthly temporal sampling. 
However, a key challenge is to discriminate 
fluctuations in ice-sheet mass from changes in the 
underlying crust and mantle. This is achieved by 
using models of GIA, which, in the case of the 
AIS, has led to significant adjustments (70). The 
spatial resolution of GRACE observations de- 
rived from global spherical harmonic solutions 
of about 300 km in the polar regions (71) is 
coarse in comparison to that of other geodetic 
techniques. Hence, a further complicating factor 
is that signals may leak into regional GRACE 
solutions as a consequence of remote geophysi- 
cal processes. In circumstances where spatial 
relationships between geophysical mass fluxes 
can be adequately characterized, application of 
the mass concentration unit (mascon) method 
(68, 72, 73) introduces a capacity to study changes 
at smaller scales. 

Data from the GRACE satellite mission were 
used to estimate changes in the mass of the AIS 
and GrIS between January 2003 and December 
2010 (16). Analysis methods varied between the 
six groups who contributed these observations; 
some used the mascon approach, whereas others 
used spatial-averaging kernels. The GRACE data 
were corrected for the effects of GIA by using the 
Simpson, ANU, and ICE-5G models in Green- 
land and the W12a, IJ05_R2, and ICE-5G mod- 
els in Antarctica. Although we only include 
results using the W12a and IJ05_R2 models in 
our reconciled estimates for Antarctica, we also 
provide separate ice-sheet mass balance estimates 
detennined by using the ICE-5G GIA model 
solution (16) to allow comparison with previous- 
ly published estimates; its use leads to more- 
negative estimates of EAIS mass balance. Each 
group made its own decisions on processing the 
GRACE data, including how to combine results 
by using different GIA models, handle contam- 
ination from external sources, compute uncer- 
tainties, and compute regional mass trends and 
time series. Ice-sheet mass time series and trends 
from all groups were then averaged to obtain the 
individual GRACE results reported here. For all 
regions, the mass trends contributed by the in- 
dividual groups agree with the combined GRACE 
trend to within the estimated uncertainties. 

Results and Discussion 

We investigated the extent to which the indepen- 
dent geodetic techniques record similar fluctua- 
tions in ice-sheet mass. First, we considered mass 
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changes within 52 AIS glacier drainage basins as 
detennined by the techniques of satellite RA and 
IOM (Fig. 1), which are well suited to this task 
( 14 , 74 ). In each case, the 19-year average rate of 
mass loss from RA was compared with values (55) 
developed by using the IOM over a similar period. 
The average difference between the estimates of 
basin mass imbalance was 1.4 + 3.8 Gtyear -1 , and 
there is agreement within 1- and 2-sigma uncer- 
tainty estimates in 42 and 49 of the 52 basins, 
respectively. Next, we computed the mass change 
as determined by satellite RA within areas of the 


EATS and WAIS that were beyond the scope of the 
IOM survey (55) to assess the extent to which 
the two methods are complementary (Fig. 1). 
These two areas, which typically fall between 
glacier drainage basins of the IOM survey, have 
small imbalances (4.5 ± 6.0 and 1 .4 + 1 .7 Gt year -1 , 
for the EATS and WAIS respectively), implying 
that the region surveyed by the IOM is sufficient 
to capture the vast majority of the present EAIS 
and WAIS mass imbalance. Furthermore, the 
IOM technique is able to resolve important mass 
changes in regions that are beyond the effective 


resolution of the RA survey, such as the APIS 
(Fig. 1). 

As a second example, we investigated the 
extent to which independent geodetic techniques 
are able to detect fluctuations in SMB. For this 
exercise, we considered an exceptional snowfall 
event in East Antarctica during the first half of 
2009 (Fig. 2). A snowfall anomaly that can be 
identified in CloudSat precipitation data (75) here 
is clearly apparent within the RACM02 (and, 
hence, IOM), RA, and GRACE data sets, which 
record the fim thickness and mass, volume, and 



Fig. 1. Comparison of ice sheet mass balance estimates derived from 
satellite RA (green) and the IOM (red) over the period 1992 to 2011, with 
1-sigma and 2-sigma error bars in dark and light shading, respectively, 
and mean values are shown in white. The comparison is performed for 


52 Antarctic drainage basins (55) and the dislocated regions of East and 
West Antarctica that are omitted from the IOM survey (EAIS_OM and 
WAIS_OM, respectively). Basin locations are illustrated in the supplementary 
materials. 





Fig. 2. Estimated anomalies in cumulative ice-sheet firn mass (A), and mass 
(B and C), derived from the RACMO regional climate model, satellite RA, and 
GRACE satellite gravimetry, respectively, over a period of anomalously high 
snowfall in East Antarctica. Anomalies were computed over the period July 
2009 to July 2010 relative to July 2008 to July 2009. Before that, linear 
trends, as fitted to the 2003 to 2008 interval, were removed. The time 
evolution of the event, as resolved by these data sets and three additional 
climate models [ERA-Interim (ERAI), CFSR, and MERRA], is also illustrated 
(D) as the average anomaly over four drainage basins of Dronning Maud 
Land in East Antarctica (shaded areas in inset map). Although there are SMB 
fluctuations elsewhere during the same time interval, the pattern of mass loss 
in West Antarctica is primarily associated with longer-term ice-dynamical im- 
balance. Relatively large annual cycles are present within some RA time 
series, but they do not obscure either short- or long-lived events, m w.e., 
meters water equivalent. 
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Fig. 3. Intercomparison of mass balance estimates of the GrIS, APIS, EAIS, WAIS, AIS, and the AIS plus 
GrIS, derived from the four independent geodetic techniques of RA (cyan), IOM (red), LA (green), and 
gravimetry (blue) over the period 2003 to 2008. Also shown is the reconciled result (gray). 




Fig. 4. Rate of mass change of the four main ice-sheet regions, as derived from the four techniques 
of satellite RA (cyan), IOM (red), LA (green), and gravimetry (blue), with uncertainty ranges (light 
shading). Rates of mass balance derived from ICESat LA data were computed as constant and time- 
varying trends in Antarctica and Greenland, respectively. The gravimetry and RA mass trends were 
computed after applying a 13-month moving average to the relative mass time series. Where temporal 
variations are resolved, there is often consistency in the interannual variability as determined by the 
independent data sets. 


mass fluctuations, respectively. The accumula- 
tion event affects the coastal region of the ice 
sheet, and, although the spatial pattern is best 
defined by the RACM02 and altimeter data sets, 
it is also apparent in the coarser-resolution 
GRACE data set. Overall, around 200 Gt of 
additional snow mass was deposited in Dronning 
Maud Land during this event, equivalent to the 
mean annual snow accumulation in this sector of 
Antarctica. In addition to the snowfall anomaly, 
the RA and GRACE data sets also include ice- 
dynamical mass changes that fluctuate over the 
survey period, such as the accelerated mass 
losses in the Amundsen Sea sector. 

In pursuit of a comprehensive methodological 
intercomparison, we computed changes in the 
mass of each ice-sheet region between October 
2003 and December 2008, the period when all 
four satellite geodetic techniques were operat- 
ing optimally (Fig. 3 and table S2). During this 
5 -year period, which is short relative to the full 
extent of the geodetic record and in comparison 
to fluctuations in SMB, the arithmetic means of 
ice-sheet mass imbalance estimates derived from 
the available geodetic techniques were -72 + 43 
and -232 ± 23 Gt year -1 for the AIS and the 
GrIS, respectively. The technique-specific esti- 
mates agree with these mean values to within 
their respective uncertainties in all four ice-sheet 
regions and for the AIS as a whole. The only 
exception is the LA estimate of the combined 
AIS and GrIS mass imbalance, which is, at 140 + 
133 Gt year -1 , more positive than the mean value 
and only marginally beyond the 1 -sigma uncer- 
tainty range of the respective values. Although 
the uncertainties of any one particular method are 
sometimes large, the combination of methods 
considerably improves the certainty of ice-sheet 
mass balance estimates. 

To produce a reconciled ice-sheet mass ba- 
lance estimate, we computed the average rate of 
mass change derived from each of the geodetic 
techniques within the various regions of interest 
and over the time periods for which geodetic 
mass rates were derived (Fig. 4). According to 
these data, ice-sheet mass balance varies cycli- 
cally and by large amounts over intermediate 
(2- to 4-year) time periods. For example, during 
the period from 1992 to 2011, the WAIS mass 
balance fluctuated around a mean value in the 
range from -50 to -100 Gt year , but there have 
been episodes of considerably larger growth and 
loss over shorter intervals. Similar variability is 
apparent in other sectors of the AIS and the GrIS. 
We next calculated the linear average of the 
individual estimates of mass balance values to 
anive at reconciled values and integrated these 
data to fonn a time series of cumulative mass 
change within each of the four ice-sheet regions 
(Fig. 5). Although there are obvious dependen- 
cies between the mass balance estimates produced 
by using each of the geodetic techniques, in- 
cluding, for example, the SMB data sets that 
are common to the IOM and LA processing, the 
G1A data sets that are common to the gravimetry 
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and altimetry processing, and the orbital cor- 
rections that are common to the LA and RA 
systems, these dependencies are in practice dif- 
ficult to characterize. For the purpose of calcu- 
lating the reconciled ice-sheet mass balance 
estimate, we considered 10M, gravimetry, and 
altimetry to be independent geodetic techniques. 
On the basis of this assumption, we compute 
the standard error of the uncertainty estimates 
from independent techniques as a measure of 
their collective uncertainty. Over the course of our 
19-year survey, the average rates of mass balance 
of the AIS and the GrlS were -71 ± 53 and -152 ± 
49 Gt year -1 , respectively (Table 1). For com- 
pleteness, we also compute cumulative mass trends 
by using the data from each individual geodetic 
technique (fig. SI). 


We also computed ice-sheet mass trends over 
shorter intervals to examine their variability (Table 
1). These estimates, along with our integrated time 
series (Fig. 5), confirm several known signals of 
mass imbalance, including increasing mass losses 
from the WAIS (55, 74 - 77 ), the APIS (75, 78 - 80 ), 
and the GrlS (5, 81 , 82 ). Although rates of mass 
loss from the GrlS were modest during the 
1990s, they have increased sharply since then 
because of episodes of glacier acceleration ( 18 , 83 ) 
and decreasing SMB (5, 66 ). GrlS glacier acceler- 
ation is, however, neither uniform nor progres- 
sive ( 65 , 84 , 85 ), and the large mass losses in 
2010 were in fact driven by anomalously low 
snow accumulation and high runoff ( 86 ). 

The WAIS has lost mass throughout the entire 
survey period, and our reconciled data set shows 


Fig. 5. Cumulative changes in 
the mass of (left axis) the EAIS, 
WAIS, and APIS (top) and GrlS 
and AIS and the combined change 
of the AIS and GrlS (bottom), 
determined from a reconciliation 
of measurements acquired by 
satellite RA, the IOM, satellite 
gravimetry, and satellite LA. Also 
shown is the equivalent global 
sea-level contribution (right axis), 
calculated assuming that 360 Gt 
of ice corresponds to 1 mm of sea- 
level rise. Temporal variations in 
the availability of the various 
satellite data sets (Fig. 4) means 
that the reconciled mass balance 
is weighted toward different tech- 
niques during certain periods. 
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Table 1. Reconciled ice-sheet mass balance estimates determined during various epochs, inclusive 
of all data present during the given dates. The period 1993 to 2003 was used in an earlier 
assessment ( 2 ). 


Region 

1992-2011 

(Gt/year) 

1992-2000 

(Gt/year) 

1993-2003 

(Gt/year) 

2000-2011 

(Gt/year) 

2005-2010 

(Gt/year) 

GrlS 

-142 ± 49 

-51 ± 65 

-83 ± 63 

-211 ± 37 

-263 ± 30 

APIS 

-20 + 14 

-8 + 17 

-12 + 17 

-29 ± 12 

-36 + 10 

EAIS 

14 ± 43 

-2 ± 54 

-9 ± 50 

26 ± 36 

58 ± 31 

WAIS 

-65 ± 26 

-38 ± 32 

-49 ± 31 

-85 + 22 

-102 ± 18 

AIS 

-71 + 53 

-48 ± 65 

-71 + 61 

-87 ± 43 

-81 ± 37 

GrlS + AIS 

-213 ± 72 

-100 ± 92 

-153 ± 88 

-298 ± 58 

-344 ± 48 


that the rate of mass loss has increased signifi- 
cantly over time (Table 1). The pattern of WAIS 
imbalance is dominated by mass losses (Amundsen 
Sea sector) and gains (Kamb Ice Stream) of dy- 
namical origin. Although close to balance during 
the 1990s, there have been significant mass 
losses from the APIS since then because of gla- 
cier acceleration in the wake of ice-shelf col- 
lapse ( 87 , 88 ) and calving-front retreat (77, 89 ). 
The APIS now accounts for around 25% of all 
mass losses from Antarctic regions that are in a state 
of negative mass balance, despite occupying just 
4% of the continental area. In contrast, the EAIS, 
which occupies over 75% of Antarctica, was in 
approximate balance throughout the 1990s. Al- 
though the EATS has experienced mass gains 
during the final years of our survey (Table 1 and 
Fig. 5), our reconciled data set is too short to 
determine whether they were caused by natural 
fluctuations that are a common feature of Antarctic 
ice-core records ( 90 ) or long-term increases in 
precipitation that are a common feature of global 
and regional climate model projections ( 91 - 93 ). 
Both satellite altimeter data sets highlight the lower 
reaches of the Cook and Totten Glaciers as re- 
gions of ice dynamical mass loss ( 15 , 77 ), but 
neither signal is large in comparison with the 
wider EAIS mass trend. Overall, snowfall-driven 
mass gains in East Antarctica, notably the 
anomalous event in Dronning Maud Land during 
2009 (Fig. 2), have reduced the rate at which 
Antarctic ice losses have increased over time, but 
the EAIS record is too short to determine whether 
this is a long-term trend. 

Our reconciliation exercise has highlighted 
several other issues. Assessments of GrlS mass 
balance require more careful consideration than 
was possible here, because the surrounding moun- 
tain glaciers and ice caps are included in some, but 
not all, of our geodetic surveys and because the 
ice-sheet domains varied in area by 2%. One esti- 
mate has put their contribution at ~20 Gt year -1 
( 94 ), a value that falls between two we have 
derived ourselves from ICESat data (10 and 40 Gt 
year -1 ). For the EAIS, our mass change estimates 
exhibit an unsatisfactory spread, with results from 
the IOM and LA techniques falling consistently 
lower and higher than the mean value we have 
derived (table S2). Although the average signal of 
EATS imbalance is relatively small, such a large 
divergence is a matter of concern; improvements 
of the ancillary data sets that support satellite ob- 
servations would be of considerable benefit in this 
region. Lastly, the spatial sampling of mass fluc- 
tuations at the APIS is at present inadequate, 
particularly considering that it provides a signifi- 
cant component of the overall AIS imbalance. 
Improvements in the spatial and temporal density 
of satellite observations of this region are needed. 

Conclusions 

We estimate that, between 1992 and 2011, the 
Antarctic and Greenland ice sheets lost 1350 ± 
1010 and 2700 ± 930 Gt of ice, respectively, 
equivalent to an increase in global mean sea level 
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of 1 1 .2 ± 3.8 mm. The greater certainty and com- 
prehensive nature of this estimate is made pos- 
sible through our combination of observations 
derived from a range of geodetic techniques, each 
of which has different strengths and weaknesses. 
We have quantified and characterized the ice- 
sheet imbalance associated with glacier dynamics 
(the APIS and WAIS), SMB (the EAIS), and a 
mixture of the two processes (the GrIS). We have 
also identified the geographical regions where 
improved data sets are required; the APIS and 
the EAIS would benefit from measurements 
with greater spatial sampling and longer tem- 
poral sampling, respectively. Although measure- 
ments from new and future satellite missions, 
such as CryoSat-2, may offer new data with which 
to tackle the former of these challenges, improve- 
ments to ancillary data sets are also required. We 
have shown that assessments of mass imbalance 
based on short geodetic records should be treated 
with care, because fluctuations in SMB can be 
large over short time periods. Lastly, our assess- 
ment demonstrates that geodetic assessments of 
ice-sheet mass balance should consider both the 
spatial and temporal limitations of the particular 
data sets upon which they are based and their 
value in relation to findings based on other inde- 
pendent approaches. 
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